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In this talk we highlight recent lattice calculations of the nucleon form factors and 
structure functions. 

t> ! 1 Introduction 

in 



The lattice formulation of QCD is, at present, the only known way of obtaining 
low energy properties of the theory in a direct way, without any model as- 
sumptions. Quantities within the grasp of lattice QCD involving light quarks 
include the hadron mass spectrum, quark masses, the A parameter, the chiral 
condensate, the sigma term, decay constants, the axial and tensor charge of 
■ the nucleon, form factors and twist two and four polarised and unpolarised 

ij moments of structure functions for the nucleon, pion and rho, to name only 

a few of the possibilities. Our group (the QCDSF collaboration) has been ac- 
| tively involved for the last few years in attempting to determine some of these 

• quantities, all characterised by their non-perturbative nature. Rather than 

giving an exhaustive progress report of recent developments in the fieldowe 
shall, due to lack of space, focus our attention on two topics: nucleon form 
factors (including the axial charge) and moments of structure functions and 
how the lattice method can lead to their determination. 

The lattice approach involves first euclideanising the QCD action and 
then discretising space-time (with lattice spacing a). The path integral then 
becomes a very high dimensional partition function, which is amenable to 
Monte Carlo methods of statistical physics. This allows correlation functions 
(which can be related to QCD matrix elements) to be determined. Progress 
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in the field is slow. First our 'box' must be large enough to fit our correla- 
tion functions into. Then, often, a chiral extrapolation must be made from a 
quark mass region around the strange quark mass to the light up and down 
quarks. Finally the continuum limit must be taken, i.e. a — > 0. This is all very 
time consuming as, in the statistical mechanics picture, we are approaching 
a second order phase transition, with all its attendant problems. Addition- 
ally, simply to save computer time, the fermion determinant in the action 
is often discarded. This 'quenched' or 'valence' procedure is an uncontrolled 
approximation. Recently, however, simulations with two mass-degenerate sea 
quarks have begun appearing, allowing a first look at the possible effects of 
quenching. Finally, in addition to the above problems, to be able to com- 
pare with phenomenological or experimental results, matrix elements must be 
renormalised. In total this is an ambitious programme. 

2 Generalities 

2.1 The processes 

Lepton-nucleon elastic scattering, IN — > IN, in which a photon is exchanged 
between the lepton (usually an electron) and the nucleon (usually a proton), 
has been studied for many years. Indeed there has been a resurgence of 
interest in these processes as part of the Jefferson Laboratory (Jlab) 'hadron' 
physics programme. The scattering matrix element can be decomposed into 
a known electromagnetic piece and an unknown QCD matrix element with 
decomposition: 

(N(p', s')\J»(q)\N(P, S)> - W, s') 

(1) 

where q = p' — p is the momentum transfer and Q 2 = —q 2 > 0. The values 
at Q 2 = are Ff (0) = 1 due to the conservation of the vector current and 
^(O) = /i p — 1, the anomalous magnetic moment, in units of e/2rriN or 
magnetons. (F™(0) = 0, ^(O) = /U™.) Experimentally, it is more convenient 
to define the Sachs form factors 

G e (Q 2 ) = F 1 (Q 2 )--^-F 2 (Q 2 ), 
(2m N y 

G m (Q 2 ) = F^Q 2 ) + F 2 (Q 2 ) . 
Earlier experimental results give phenomenological (dipole) fits of 

Gl{Q 2 ) = = = ( X + 5) 2 ' G " (Q2) = ' (2) 
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7 M ^i(Q 2 )+^ 



2m n 



F2(Q 2 ) 



u{p, s) , 



with m v w 0.82 GeV, /J," w 2.79 and //> w -1.91. 

Similarly neutrino-nucleon scattering, for example v^n — > mediated 
by a H^" 1 " exchange, leads to an unknown axial current hadronic matrix element 
between neutron and proton states, which, with the use of current algebra and 
isospin invariance, may be re- written between p states alone and has Lorentz 
decomposition 

(pip', s")\A^ d (q)\p(p, 5)) = u(p\ s') 

(3) 

where A™ d — uj^j^u — drf^sd. Experimental results give phenomenological 
fits 

9 a(Q 2 )=9a(0)(i + ^-) , (4) 

V m A/ 

with m A w 1.00 GeV. From the (3 decay we know! g A = g A (0) = 1.2670(35). 

At higher momentum transfer, the nucleon is broken up by the photon 
(or W ) probe. We enter the regime of deep inelastic scattering experiments 
eiV — > eX (or v^n — > li~X). The operator product expansion, OPE, leads 
to relations between moments of the structure functions and certain nucleon 
matrix elements. For example, 

j\xx n - 2 F 2 {x,Q 2 ) = l E ^(^,g^)v^(gn + 0(^). (5) 

Here x is the Bjorken variable and v" s oc (N\O n \N) MS , where O n ~ 
g7 Ml £>^ 2 . . . D fin q are operators bilinear in the quarks, each containing n — 1 
covariant derivatives. For e./V scattering only even moments contribute, while 
for vN all moments are allowed. 

All the matrix elements briefly discussed above can be determined non- 
perturbatively using lattice QCD. 



u(p,s) 



2.2 Lattice Technicalities 

The general method of determining matrix elements is by forming ratios of 
three-point to two-point correlation functions: 

^ ' (N(t;p)N(0;p)) K ' 

where N is the three quark baryon operator and O is a bilinear quark operator. 
Provided that f » t > 0, then R can be shown to be directly proportional 
to the matrix element (N\0\N). (R may be easily generalised to allow for 
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momentum transfer.) If we re- write the numerator of eq. (0) using quark 
propagators, we see that we have two classes of diagrams: one in which the 
inserted operator is attached to quark propagators to the nucleon, and a 
second class where the quark propagator from the operator is disjoint from 
the baryon. In this second class the interactions between the inserted operator 
and the baryon take place only via gluon exchange. Numerically, due to large 
ultra-violet fluctuations, no useful signal is seen. So either we must presently 
assume that their effects are small, or calculate 'non-singlet' (i.e. u — d) matrix 
elements where these second class terms cancel. 

The matrix elements must be renormalised. There are several possibil- 
ities. We may simply consider a ratio of matrix elements between different 
hadronic states, so that the renormalisation constant cancels. However, this 
is not very useful. The axial and vector renormalisation constants (the cur- 
rents being partially conserved) may be determined by demanding that their 
(continuum) Ward Identities are obeyed, which produces non-perturbative 
renormalisation constants. (Incidently they are also renormalisation scheme 
independent.) Perturbation theory can be applied, but due to technical prob- 
lems only one loop results are known. Even with this restriction, perturbation 
theory can be improved leading to 'tadpole improved' or TI perturbation the- 
ory. However it still suffers from unknown systematic errors. The ALPHA 
Collaboration has developed techniquestJ allowing some renormalisation con- 
stants to be determined by stepping up from a low energy scale to a high 
energy scale where the Z can be matched to (known) perturbation theory. 
Finally one can attempt to 'mimic' the perturbation theory method by cal- 
culating the operator numerically between quark states in the Landau gauge. 
Also note that it is desirable to make the errors in the determination to be 
of 0(a 2 ) rather than 0(a), because then the approach to the continuum is 
fasteru. This introduces further irrelevant operators, with 'improvement' co- 
efficients which also must be determined. 

At present, for the axial and vector currents most renormalisation con- 
stants and improvement coefficients are known non-perturbatively. For v n we 
rely on TI perturbation theory. 

3 Results 

3.1 The Axial Charge 

Using the quenched approximation we have made simulations at three lattice 
spaemgs: a = 0.093, 0.068 and 0.051 fmQ. For each lattice size simulations are 

a Corresponding to a" 1 = 2.12,2.90 and 3.85 GeV, with scale! r = 0.5fm. 
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made at four or more quark masses, which are then linearly extrapolated to the 
chiral limit. (A picture for v n will be shown later in section 3.3.) We also haw; 
one unquenched result at a rs 0.11 fm. In Fig. |l| we show all these numbersu. 
The results all lie somewhat lower than the experimental value, and at present 
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Figure 1. The continuum extrapolation of gA- 



there is no clear difference between quenched and unquenched fermions. (Note 
that for the quenched results, this is a fully non-perturbative calculation; for 
the unquenched result, there is still a little uncertainty.) While it is possible 
to extrapolate the quenched data to the continuum limit (although it would 
be very desirable to have more points because the noise is uncomfortably large 
at present) , we see that we are only at the beginning of obtaining results with 
dynamical fermions. 

3.2 Form Factors 

In Fig. |^ we show proton electric and magnetic form factors^, for our quenched 
lattice a — 0.068fm. Other lattice values and the unquenched results are 
very similar. Due to the momentum transfer the results are rather noisy. 
However the present trend is clear: the dipole fit gives values of my too 
large in comparison with the phenomenological value. This is confirmed in 
Fig. ^. For the nucleon we find r rms = y/l2/mv ~ 0.60 fm in the quenched 
approximation and r rms w 0.70 fm for the dynamical case. This is to be 
compared with the phenomenological value r rms = 0.83 fm. Perhaps the 'pion 
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Figure 2. and Gj, compared to experimental data for quenched fermions at a = 0.068fm. 
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Figure 3. m„ e compared to its phcnomcnological value. 



cloud' has not developed at our dynamical quark masses, which makes the 
nucleon appear to be smaller than it really is. pP is roughly consistent with 
the phenomenological value (but with large errors), while trip™ is again too 
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large. _ 

The recent Jlab resultsa for Gf/G^ indicate a functional difference be- 
tween the electric and magnetic form factors. This is shown in Fig. H. If 
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Figure 4. /GjJ, compared to experimental data for quenched fermions at a = 0.068fm. 

cq. (||) holds then the ratio of form factors should be constant. Unfortunately 
at present the lattice data is too noisy to draw any conclusion. 

Finally in Fig. ^| we show the axial form factor. Again the previous com- 
ments as for Gf , G v m hold verbatim. 

3.3 Moments of the nucleon structure function 

Some years ago, we published our first results! for v n , n = 1, 2, 3. The re- 
plotted data for v^ MS — Vn is shown again in Fig. ^. for unimproved Wil- 
son fermions. A comparison is made with MRS phenomenological resultscEl. 
While higher moments tend to agree (although large error bars might hide any 
discrepancy), it seems that there is a discrepancy for the lowest moment. This 
may be due to renormalisation and discretisation effects, higher- twist contri- 
butions, the chiral extrapolation^, or perhaps the phenomenological distri- 
bution functions, being global fits, do not fit this moment very well. Naively 
this result is also perhaps expected as V2 is part of the energy-momentum 
sum rule, ^■ i; ™ lence _|_ ^| ea = i ; anc [ d ue to quenching the sea contribution is 
reduced. For higher moments this is less of an effect, as the gluon contribution 
is more important for smaller x. 
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Figure 5. g P A compared to experimental data for quenched fermions at a = 0.068fm. 
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Figure 6. The chiral extrapolation for v^ MS — v!$^i S for quenched fermions at a = 
0.093fm (circles), compared to the MRS fit functiontJ (filled squares). The dot-dashed 
lines show the positions of (% and a (hypothetical) strange quark pseudoscalar particle. 



First let us look at the experimental data. The cleanest direct determina- 
tion of v^ Mb — vffl MS is given from the F$ — F£ NMC dataS, as shown in 
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Figure 7. Thf-experimental binned data for F% - at Q 2 ] = 4GeV 2 from the NMC 
Collaborationt3 compared with the MRS parton fit function . 



Fig. [7], because the one experiment uses both p and n targets. Note that there 
is a paucity of data at large x ~ 1, which is the region where the moments are 
most sensitive to. (We make a simple linear extrapolation from the largest 
data point to x = 1.) While the MRS fit does not quite follow the data, the 
area under the curve cx v^ MS — v^" 3 is about the same, so there is indeed 
no problem. 

Secondly we have investigated the continuum extrapolation for the lowest 
moment, in Fig. |^ using improved fermions. Surprisingly, perhaps, there seems 
to be little a 2 trend in the data (c.f. this with g^, Fig. |l|). 

Thus the previous discrepancy still holds. 
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